Universal optimal transmission of light through disordered materials 
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We experimentally demonstrate increased diffuse transmission ol light through strongly scattering 
materials. Wavefront shaping is used to selectively couple light to the open transport eigenchannels, 
specific solutions of Maxwell's equations which the sample transmits fully, resulting in an increase 
of up to 44% in the total angle integrated transmission compared to the case where plane waves 
are incident. The results for each of several hundreds of experimental runs are in excellent quan- 
titative agreement with random matrix theory. From our measurements we conclude that with 
perfectly shaped wavefronts the transmission of a disordered sample tends to a universal value of 
2/3, regardless of the thickness. 

PACS numbers: 05.60.-k, 42.25.Dd, 73.23.-b 
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The transport of waves in strongly scattering media is, 
usually, well described by diffusion theory. However, the 
diffusion equation does not take into account interfer- 
ence. Interference gives rise to fundamental effects such 
as enhanced backscattering[l|, Q and Anderson localiza- 
tion of light HH i|. Such interference effects may be ob- 
served with incident light in a single freely propagating 
mode (scattering channel), such as a plane incident wave. 
In 1984, Dorokhov predicted a striking multichannel in- 
terference effect. Using random matrix theory (RMT), he 
showed that the transmission through a diffusive mate- 
rial is the result of a small number of open eigenchannels 
with a transmission coefficient of close to one. [6] Each 
eigenchannel [l!| corresponds to a specific linear combi- 
nation of multiple free modes. Waves coupled to an open 
eigenchannel will be fully transmitted through a disor- 
dered sample, even if the sample is optically thick. In 
contrast, an incident plane wave couples mainly to closed 
eigenchannels (with a transmission coefficient of close to 
zero) and, consequently, most of its power is diffusely re- 
flected. This result was originally obtained for electron 
transport in a wire. Later, it was generalized to a slab 
geometry and to optical systems [1, [§]■ 

Universal conductance fluctuations, which are ran- 
dom fluctuations in the coupling of waves to open and 
closed eigenchannels, have been observed in experiments 
13, EI El El, however, no controlled coupling to open 
eigenchannels has been reported so far. Here we experi- 
mentally demonstrate the injection of light into the open 
eigenchannels of a strongly disordered sample, resulting 
in a large increase of the diffusely transmitted intensity. 
We accomplish this by shaping the wavefront of the inci- 
dent light, which amounts to individually controlling the 
phases of multiple incoming free modes. 

The transmission amplitudes of waves through a meso- 
scopic medium are described by a transmission matrix t. 
The field in the outgoing free modes at the back of the 




FIG. 1: Experimental setup. HeNe, expanded 632.8 nm HeNe 
laser; A/2, half waveplate; PBS, polarizing beam splitter cube; 
SLM, spatial light modulator. A, iris diaphragm; 63x, micro- 
scope objective; 60x, oil immersion microscope objective; S, 
sample; P, polarizer; ND, neutral density filter; Li, L2, L3, 
lenses with focal length of respectively 250 mm, 200 mm and 
600 mm. ND and S are translated by computer controlled 
stages. 



sample is given by 
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where indices a and b label incident and transmitted free 
modes respectively. N is the total number of incident 
free modes. Any complete orthogonal set of modes may 
be chosen; in mesoscopic physics it is usual to choose 
the transversal modes of a perfect waveguide as a basis, 
one may also work with an overcomplete basis of incident 
angles or diffraction limited spots on the sample surface. 
Transmission eigenchannels are defined for each individ- 
ual sample by decomposing the transmission matrix as 
t = UTV T . Here U and V are unitary matrices that ef- 
fect the basis transformations between free modes outside 
the sample and eigenchannels inside. T is a real diagonal 
matrix containing eigenchannel transmission coefficients. 

We use a spatial light modulator (see Fig. Q]) to shape 
the wavefront of the incident light. The computer pro- 
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gram that controls the wavefront shaper optimizes the 
intensity of a small target spot at the back surface of 
the sample using the measured intensity in the target as 
feedback This wavefront shaping method creates a 
high intensity focus through the sample, which by itself is 
not a result of the bimodal distribution of the eigenchan- 
nel transmission coefficients. Importantly, the disordered 
background speckle around the focus turns out to be a 
very sensitive probe of the distribution of the elements of 
T. 

For example, for the limiting distribution were all 
eigenchannels are completely open T = or completely 
closed T = 1 (the so called maximal fluctuations [H), 
only open eigenchannels will contribute to the intensity 
in the target. Therefore, the optimization algorithm will 
select only linear combinations corresponding to open 
eigenchannels. As a result, the diffuse transmission will 
increase so that the ideally shaped wavefront will have 
a total transmission of unity. The distribution predicted 
by DorokhovQ is close to this limiting case. 

We now proceed to introduce a quantitative measure 
of the control we exert over a shaped wavefront. Our 
algorithm maximizes the intensity in a diffraction limited 
spot, which is exactly one of the transmitted free modes. 
We label this special target mode with the index (3. The 
ideally shaped incident wavefront E° pt '@ for maximizing 
the intensity in /3 is given by [lij], 
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where Tp normalizes the total incident power. Our op- 
timization algorithm proceeds as follows: the matrix el- 
ements tp a are measured up to a constant prefactor by 
cycling the phase of the light in the incident mode a while 
observing the intensity in target mode (3 14J. After N 
phases have been measured, the optimized incident wave- 
front is constructed according to Eq.Q. This optimized 
wavefront couples to a superposition of eigenchannels, 
mostly to channels with high transmission eigenvalues. 

In any experiment, the resolution and the spatial ex- 
tent of the generated field are finite. Therefore, it will 
never be possible to exactly construct the wavefront de- 
scribed by Eq. |2]). To quantify how well the actual 
incident field E^ ct '^ matches the optimal incident field 
j^opt,/^ we introduce the overlap coefficient 7 as 
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The degree of intensity control is I7I 2 . We can now write 
any incident wavefront as a linear superposition of the 
perfect wavefront and an error term 

^iact,/3 
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(4) 



where the error term AE a is normalized. For ideal con- 
trol over the incident wavefront I7I 2 = 1. In earlier 



experiments 14| the degree of control was relatively low 
(I7I 2 <C 0.1), and total transmission did not increase mea- 
surably. In this Letter, we discuss experiments at much 
higher values of I7I 2 , up to 0.33. 

The experimental apparatus (see Fig. [T]) is designed 
to approach the optimal wavefront as closely as possible, 
by controlling the largest possible fraction of the inci- 
dent free modes. An expanded beam from a 632.8 nm 
HeNe laser is rotated to a 45° linear polarization by a 
half waveplate and impinges on a polarizing beam split- 
ter cube. Horizontally and vertically polarized beams are 
modulated with separate reflective liquid crystal displays 
(Holoeye LC-R 2500) and then recombined, to provide 
control over modes with both polarizations. We used a 
4- pixel macropixel modulation method [lj| to control the 
phase of the light without residual amplitude modula- 
tion. The modulator is divided into 3816 independently 
programmable segments. A computer programs the mod- 
ulators using feedback from a camera as discussed below. 
A sequential optimization algorithm 14] was used to op- 
timize the wavefront [l5} . A high numerical aperture ob- 
jective (NA = 0.95, Zeiss Achroplan 63 x) projects the 
shaped wavefront onto the sample. 

Each sample consists of a layer of spray-painted ZnO 
particles on a standard glass microscope cover slip. The 
particles have an average diameter of 200 nm, which 
makes them strongly scattering for visible light. The 
mean free path was determined by measuring the total 
transmission and equals 0.85 ± 0.15 /jm at a wavelength 
of 632.8 nm. We used samples with thicknesses of 5.7 /j,m 
and 11.3 /im. The samples were positioned in the focal 
plane of the microscope objective to minimize the size of 
the diffuse spot, and thereby the number of contributing 
modes. The number of such contributing free modes was 
estimated from the intensity profile of the transmitted 
light [lH to be 5.5 • 10 3 and 1.0 • 10 4 modes, for the thin 
and the thick samples, respectively. The samples were 
mounted on a motorized stage to translate them in the 
focal plane. 

A high NA oil-immersion objective (Nikon TIRF 
60x/NA= 1.49) collects the transmitted light. The 
transmitted light is split into horizontal and vertical po- 
larizations by a beam splitter cube. A second polarizer 
improves the extinction ratio for reflected light. The 
magnification of the detection system is 225 x , enough to 
well resolve individual speckles. A camera measures the 
power of the horizontally polarized light in a disc with a 
diameter of 0.11 fim. at the sample, which is smaller than 
a single speckle, to provide feedback for the optimization 
algorithm. 

After optimization, a calibrated neutral density filter 
with a transmission of 1.4 • 10~ 3 is placed in front of 
the camera to measure the high intensity in the target. 
A second camera images the intensity of the vertically 
polarized light. 

Optimizing the incident wavefront caused the intensity 
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FIG. 2: Intensity distribution of horizontally polarized light 
in a 15 /»m x 15 /im square area at the back of the sample, (a) 
For a non-optimized incident wavefront. (b) For an optimized 
wavefront. (c) Intensity summed in the y-direction to average 
over speckle. Dashed curve, transmission of non-optimized 
wavefront; Solid curve, transmission of optimized wavefront. 



in the target to increase dramatically. In Fig. [2] we plot 
the transmitted intensity through a 11.3 /zm-thick sam- 
ple for a non-optimized wavefront and for the optimized 
wavefront 15j. Before optimization, the transmitted in- 
tensity forms a diffuse spot on the back surface of the 
sample. After optimization, a strong peak emerges in 
the target focus. The intensity increase in the center of 
the target was a factor of 746 ± 28. After optimization, 
2.3% of the incident power is transmitted into the target 
focus. 

More importantly, the intensity in an area with a ra- 
dius of approximately 5 /im around the target also in- 
creased, even though the algorithm did not use this in- 
tensity as feedback. This observation indicates that we 
have redistributed the incident light from closed eigen- 
channels to open eigenchannels. As a result of optimizing 
a single target point, the total angle-integrated transmis- 
sion increased from 0.23 to 0.31. This change amounts 
to a relative increase of 35%. 

For a quantitative analysis we need to know the degree 
of control \"f\ 2 ■ Factors like measurement noise and ther- 
mal drift result in a different degree of control for each 
single run of the experiment. Fortunately, it is possible 
to measure 7I 2 directly for each run by observing the 
intensity in the target. Only the controlled fraction of 
the incident wavefront contributes to the intensity in the 
target. The transmission to the target mode (3 equals 
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where we used Eq. (jj) and the fact that the error term 



is orthogonal to the ideal wavefront. 
Eq. ([2]) we obtain 



By substituting 
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\Y\Tp- (7) 
Equation ([7]) allows us to obtain the degree of control 

2 

by measuring the intensity in the target focus Ep ct 

and Tp. In the experimental procedure it is very im- 
practical to measure Tp. Therefore, we approximate 
Tp = Tt tN/M. Here, T to t is the ensemble averaged total 
transmission of an unoptimized wavefront, and M is the 
number of transmitted free modes. Since our samples are 
sandwiched between a glass substrate on one side and air 
on the other side, the number of modes on the back of 
the sample is larger and M = n 2 N, with n = 1.52 the 
refractive index of the substrate. This approximation ne- 
glects the Ci fluctuations in the total transmission, 
which are in the order of 2% for our samples. 

For the experimental run that is shown in Fig. [2] we 
find a degree of control of |7| 2 = 0.23. This means that 
the incident field is a linear superposition of the perfectly 
shaped wavefront (carrying 23% of the incident power) 
and a noise term (carrying the rest of the power). The 
total, angle-integrated transmission T^ contains contri- 
butions both from the perfectly shaped wavefront and 
from the noise term, 

T£ = T c + (1 - | 7 | 2 ) Tt„t, (8) 

where T c is the part of the transmission resulting from 
the perfectly shaped fraction of the incident wavefront. 
By substituting Eq. @ into Eq. (fTJ) and summing the 
power in all transmitted free modes, we find 
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We evaluate 64,2 theoretically by averaging over all pos- 
sible target modes j3. We assume that C4 2 is self aver- 
aging, which is verified by our experiment. Neglecting 
small correlation terms between numerator and denomi- 
nator we find 
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From Eq. (TIT))) , it becomes clear that C4.2 is a mea- 
sure for the width of the distribution of the transmission 
eigenvalues By measuring the total transmission after 
optimizing the incident wavefront, we have direct exper- 
imental access to this value for each single sample. Since 
we measured \j\ 2 separately, we can use Eqs. (JSJ and © 
to obtain 64,2 from a single, non-ideal experimental run. 
In the particular run in Fig. [2l we find (7^2 = 0.62. 

The ensemble averaged value for C^^ was derived us- 
ing RMT. RMT[3| predicts (C 4i2 ) = 2/3 for a non- 
absorbing system far away from the localization transi- 
tion, regardless of the original transmission coefficient of 
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FIG. 3: Total transmission of the controlled fraction of the 
wavefront. Data from two different samples collapses to a 
universal curve. There are no adjustable parameters in the 
data processing or the theory. Circles, results for 11.3 /xm- 
thick sample; squares, results for 5.7/im-thick sample; star, 
data point corresponding to the experimental run that was 
shown in detail in Fig. [2] solid line, linear regression with a 
slope of 0.68; dotted line, slope of 2/3 as predicted by RMT. 



the system. For a single realization of disorder, we found 
Ci,2 = 0.62. To investigate the universality of this re- 
sult and to compare the measured values with RMT, we 
performed automated sequences of measurements. The 
sample was translated between the runs to obtain differ- 
ent realizations of disorder. 

The intensity increase was different for each of the opti- 
mizations in the sequence. The increase in the integrated 
transmission varied from a few percent to a maximum of 
44%. These variations are the result of drift due to vary- 
ing environmental conditions. While drift is undesirable, 
it gives us a wide range of |7| 2 to investigate. The degree 
of control I7I 2 was determined for each of the runs using 
Eq. 0. Then, using Eq. ©, we isolated T c . 

In Fig. [3] we plotted T c versus \j\ 2 for hundreds of 
measurements. To determine the effect of the sample 
thickness, two different samples were used: one with a 
thickness of 5.7 ± 0.5 /im, and one that is approximately 
twice as thick (11.3 ± 0.5 /urn). All data points collapse 
to a single line. Our data shows that the transmission 
coefficient of an ideally shaped wavefront does not de- 
pend on the sample thickness. A linear regression gives 
(C4,2) = 0.68 ± 0.07 where the uncertainty follows from 
a worst case estimate of the systematical errors in the 
experiment. This value is in excellent agreement with 
RMT. 

We have shown experimental evidence of controlled 
coupling of light into open transmission eigenchannels in 
opaque, strongly scattering materials. The coupling to 
open eigenchannels was enhanced by using a wavefront 
shaping algorithm to optimize transmission to a focus, 



and detected by measuring the angle-averaged transmis- 
sion intensity, which showed a relative increase of up to 
44%. We quantitatively compared the results for different 
samples and different realizations of disorder. All results 
showed a universal behavior that is in excellent quantita- 
tive agreement with random matrix theory. Our results 
demonstrate that RMT of wave transport can success- 
fully be applied to open systems and single realizations 
of disorder. This conclusion is relevant for the propaga- 
tion of electromagnetic waves, matter waves and sound 
in open, strongly scattering environments. 
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